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ISOTROPIC 7-ray BACKGROUND: COSMIC-RAY INDUCED ALBEDO FROM DEBRIS IN THE SOLAR SYSTEM? 
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ABSTRACT 

We calculate the 7-ray albedo due to cosmic -ray interactions with debris (small rocks, dust, and grains) in the 
Oort Cloud. We show that under reasonable assumptions a significant proportion of what is called the "extra- 
galactic 7-ray background" could be produced at the outer frontier of the solar system and may be detectable 
by the Large Area Telescope, the primary instrument on the Fermi Gamma-ray Space Telescope. If detected 
it could provide unique direct information about the total column density of material in the Oort Cloud that is 
difficult to access by any other method. The same 7-ray production process takes place in other populations of 
small solar system bodies such as Main Belt asteroids, Jovian and Neptunian Trojans, and Kuiper Belt objects. 
Their detection can be used to constrain the total mass of debris in these systems. 

Subject headings: interplanetary medium — Kuiper Belt — Oort Cloud — cosmic rays — diffuse radiation 
— gamma-rays: theory 



1. INTRODUCTION 

Studies of the outermost region of the solar system started 
in antiquity as documented by the historical record of comet 
observations (Kronk 1999). Among the oldest are Babylo- 
nian inscriptions referring to a comet of 674 B.C. The first 
proven observation of the IP/Halley comet is recorded in the 
Chinese text Shih clii dated 239 B.C., and its revisits of the 
inner solar system have been recorded by subsequent genera- 
tions of astronomers. However, ancient astronomers had very 
little knowledge about the origin of these comets. Only in 
1950, based on observations of a handful of so-called long- 
period (LP) comets was the conclusion made that there is a 
huge reservoir of ^10^^ comets that surrounds the solar sys- 
tem at distances larger than '^10'* AU (Oort 1950). Since that 
time the total number of identified LP comets has increased 
to ~400 (Marsden & Williams 2003). This reservoir is now 
known as the Oort cloud (OC) and its outer edge is placed be- 
tween 5 X 10"* and 2 x 10^ AU 1 pc (Dones et al. 2004; 
Weissman 2000). The OC population is the least explored in 
the solar system because its extremely distant location from 
the Sun renders the detection methods used for objects such 
as Kuiper Belt Objects (KBOs) ineffective^. 

Perturbations due to the Galactic tide, passing stars, and gi- 
ant molecular clouds alter the orbits of objects in the OC. The 
result is that a small fraction of these objects are injected into 
the inner solar system, appearing as LP comets. Such comets 
are the only members of the OC population that are available 
for study by astrophysical methods so far. Therefore, conclu- 
sions on the structure of the OC and composition of its bodies 
are mainly based on the studies of LP comets. As the comets 
are largely composed of ices, other OC bodies are assumed 
to have a similar composition. The sizes of the comet nuclei 
range between a few km to ^50 km. A population of bodies 
below a sub-kilometer size is terra incognita. 

The OC is believed to be a remnant of the proto-planetary 

' Also Kavli Institute for Particle Astrophysics and Cosmology, Stanford 
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^ Here we mean the usual optical- and infrared-based methods relying on 
the detection of scattered solai' radiation from the individual objects. 



disk that formed around the Sun approximately 4.6 billion 
years ago. The total mass and geometry of the OC are deter- 
mined by very indirect methods and rely mostly on computer- 
intensive simulations (e.g., Dones et al. 2004). The planetesi- 
mals comprising the OC initially coalesced much closer to the 
Sun, but gravitational interactions with giant planets ejected 
the proto-planetary material into the outer edge of the solar 
system (e.g., Morbidelli 2005). As the accretion and collision 
rates in the OC are extremely slow, it is believed that it con- 
tains the most pristine material left over from the epoch of 
planet formation. 

A new method to study the populations of small solar sys- 
tem bodies (SSSBs - such as Main Belt asteroids (MBAs), 
Jovian and Neptunian Trojans, and KBOs) has been recently 
proposed by Moskalenko et al. (2008). It was shown that 
the cosmic -ray (CR) induced 7-ray flux from SSSBs which 
are large enough for the CR cascade to fully develop in the 
rock (typically > 1 meter in diameter - the "thick target" 
case) strongly depends on the size distribution of the SSSB 
populations and may be detectable with the Large Area Tele- 
scope (LAT), the primary instrument on the Fermi Gamma- 
ray Space Telescope (Fermi) (formerly the Gamma-ray Large 
Area Space Telescope (GLAST)). The 7-ray spectrum from 
these processes is very steep with an effective cutoff around 
1 GeV. If detected by the LAT, it can provide unique infor- 
mation about the number of objects with sizes down to a few 
meters in each system. 

For objects with diameters below ^ 1 meter, the CR cas- 
cade does not fully develop since the total column density of 
material is < 1 interaction length. In this "thin target" case, 
the albedo 7-rays depend only on the total column density of 
material in a particular direction. The albedo spectrum has a 
shape characteristic of pp-interactions where the high energy 
spectral slope is similar to that of the incident CR spectrum. 
Therefore, it is possible also using 7-ray observations to put 
constraints on the total amount of debris in each system. 

The outermost frontier of the solar system, the OC, has a 
wide distribution on the sky. This is different to other popula- 
tions of small bodies such as the MBAs, Jovian and Neptunian 
Trojans, and KBOs which are distributed near the ecliptic. 
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The albedo 7-rays produced by CR interactions with debris 
in the OC apparently contribute to the isotropic 7-ray back- 
ground (IGRB), which is usually assumed to be extragalactic. 
A correct estimate of the 7-ray albedo of the OC is important 
for disentangUng the truly extragalactic component. On the 
other hand, an estimate of the total column density of the dis- 
persed material in the OC may shed light on the early history 
of the solar system and constrain planetary evolution models. 

2. ISOTROPIC 7-RAY BACKGROUND 

The IGRB was first discovered by the SAS-2 satellite 
(Thompson & Fichtel 1982) and confirmed by EGRET 
(Sreekumar et al. 1998). The IGRB is thought to be a su- 
perposition of all unresolved sources of high-energy 7-rays in 
the universe plus any truly diffuse component. A list of the 
contributors to the IGRB includes "guaranteed" sources such 
as blazars and normal galaxies (Bignami et al. 1979; Pavlidou 
& Fields 2002), and potential sources such as galaxy clusters 
(Ensslin et al. 1997), shock waves associated with large scale 
cosmological structure formation (Loeb & Waxman 2000; 
Miniati 2002), distant 7-ray burst events (Casanova et al. 
2007), pair cascades from TeV 7-ray sources and ultra-high 
energy CRs at high redshifts (so-called Greisen-Zatsepin- 
Kuzmin cut-off). A consensus exists that a population of un- 
resolved AGN certainly contribute to the IGRB inferred from 
EGRET observations; however predictions range from 25% 
up to 100% of the IGRB (Chiang & Mukherjee 1998; Mucke 
& Pohl 2000; Mukherjee & Chiang 1999; Stecker & Sala- 
mon 1996). A summary of the conventional contributors to 
the IGRB can be found in Dermer (2007). A number of ex- 
otic sources that may contribute to the IGRB have also been 
proposed: baryon-antibaryon annihilation phase after the Big 
Bang, evaporation of primordial black holes, annihilation of 
so-called weakly interacting massive particles (WIMPs), and 
strings. 

The IGRB is a weak component which is difficult to disen- 
tangle from the intense Galactic foreground. Extensive work 
has been done (Sreekumar et al. 1998) to derive the spec- 
trum of the IGRB based on EGRET data. A new detailed 
model of the Galactic diffuse emission (Moskalenko & Strong 
2000; Strong et al. 2004a) lead to a new estimate of the IGRB 
(Strong et al. 2004b) which is lower and steeper than found 
by Sreekumar et al. (1998); it is not consistent with a power- 
law and shows some positive curvature, as expected for an 
origin in blazars. But, recent work has shown deficiencies in 
the understanding of the IGRB. Two more diffuse emission 
components originating nearby in the solar system have been 
identified: 7-ray emission due to inverse Compton scattering 
of solar photons by CR electrons (Moskalenko et al. 2006; 
Orlando & Strong 2007) and a 7-ray glow around the ecliptic 
due to the CR-induced 7-ray albedo of SSSBs (Moskalenko et 
al. 2008). The former has been detected in the EGRET data at 
a level consistent with predictions (Orlando & Strong 2008). 
Significantly, the level of emission from this process is com- 
parable to the current estimate of the IGRB within 10° of the 
Sun path on the sky. 

3. THE OORT CLOUD 

Although an extensive literature on the origin, population, 
and dynamics of the OC exists, properties of the OC popula- 
tion are derived mostly from computer simulations and stud- 
ies of the LP and Halley-type comets which are thought to 
originate in the OC. In this section, we give a brief overview 
of the structure of the OC relevant to the current investigation. 
For more detailed information we refer the reader to the ex- 
haustive books and reviews (e.g., Brandt & Chapman 2004; 



Dones et al. 2004; Fernandez 2005, and references therein). 

Observation of a sharp spike in the number distribution of 
comets at near-zero but bound energies, representing orbits 
with semimajor axes exceeding 10^ AU, led Oort to conclude 
that the spike, a huge near- spherical cloud of icy bodies at 
> 2 X 10* AU (outer OC), had to be the source of the LP 
comets. Hills (1981) has shown that the apparent inner edge 
of the OC at a semimajor axis a ^ (1 — 2) x 10^ AU could be 
a selection effect due to the small probability of close stellar 
passages capable of perturbing comets at smaller distances. 
The number of comets that reside in the Hills cloud (inner 
OC) at semimajor axes of a few thousand AU could be signif- 
icantly larger than in the outer OC. The inner OC could also 
serve as a reservoir that replenishes the outer OC stripped by 
an external perturber Close passages of stars perturbing the 
inner OC could result in so-called "comet showers", the ex- 
treme increases in cometary flux which lasts for a few orbital 
periods after the passage. One such shower may be respon- 
sible for a ^^2.5 My period of increased bombardment of the 
Earth which produced the large Popigai (100 km) and Chesa- 
peake Bay (90 km) craters and several smaller craters ^36 
million years ago (Farley at al. 1998). 

A possible clue to the existence and the structure of the in- 
ner OC could be the Halley-type comets. Dynamical sim- 
ulations show that the source of these comets (Levison et 
al. 2001) is required to be a massive doughnut- shaped inner 
OC with median inclination between 10° and 50°. Yet, there 
could exist a third innermost region, the so-called inner core 
(Fernandez 2005), which may span between the Scattered 
Disk component of the Kuiper Belt and the inner OC. The 
inner core may harbor a large number of comets and mass, 
but does not exhibit itself due to the lack of very close stel- 
lar passages. Dwarf planets 90377 Sedna and 148209 2000 
CR105 could be the first inner core candidates. 

Dynamical models put the total number of comets with a > 
(1 — 2) X 10^ AU at ~10^^. Assuming an average mass for 
a comet 4 x 10^^ g, the total mass is estimated at 7Mq . The 
more massive inner OC could harbor ^^(2 — 13) x 10^^ comets, 
which yields a mass of 14 — 90 M©. There is no information 
on the total mass and number of sub-kilometer size bodies. 

Even closer to the Sun, beyond Neptune's orbit there is the 
Kuiper Belt. The KBOs are not uniformly distributed, with 
at least three dynamically distinct populations identified: the 
Classical Disk, the Scattered Disk with large eccentricities 
and inclinations, and "Plutinos" around the 3:2 mean motion 
resonance with Neptune at 39.4 AU. KBOs are distributed 
between 30 - 100 AU (Backman et al. 1995, and references 
therein). A majority of the KBOs have their orbits distributed 
near the ecliptic with FWHM of the order of 10° in ecliptic 
latitude (Brown 2001). The total mass of KBOs is estimated 
to be in the range ~0.01-0.3 M©, while the most often used 
value is ^0.1 Mq (Luu & Jewitt 2002). A summary of the 
mass and size distributions of KBOs, MBAs, and population 
of small bodies in Jovian and Neptunian Trojan famihes can 
be found in Moskalenko et al. (2008). 

The mass and size distributions for populations of asteroid 
families or other SSSBs are thought to be governed by colU- 
sional evolution and accretion. Under the assumptions of scal- 
ing of the collisional response parameters and an upper cutoff 
in mass, the relaxed size distributions approach a power-law 
dN OC r~"dr (Dohnanyi 1969), where r is the radius of the 
body, and n = 3.5 for a pure Dohnanyi cascade. In reality 
different populations of SSSBs show deviations from this in- 
dex. 

The dust particles and grains in the solar system are sub- 
ject to many processes which produce a negligible effect on 
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larger bodies, therefore simulations of the dynamics of the 
asteroid populations can not tell much about the total mass 
of debris nor the size distribution. The main forces acting 
on small particles are (e.g., Griin 2000): gravitation, radi- 
ation pressure which is directed outward from the Sun, the 
Poynting-Robertson effect which is essentially a drag reduc- 
ing the orbital eccentricity and causing the particles to slowly 
spiral toward the Sun, sublimation, collisions among dust par- 
ticles, plus interactions with the solar wind and magnetic field 
due to a small electric charge carried by grains which fur- 
ther complicates their dynamics. The relative effect of these 
processes is different across the range of particle sizes and 
number densities. Thus, it is very difficult to simulate the dis- 
tribution of debris in the solar system. However, if the column 
density of the debris is large enough its distribution can be di- 
rectly probed with the FermZ/LAT 7-ray telescope. 

4. CALCULATIONS 

Our previous calculations (Moskalenko et al. 2008) consid- 
ered bodies large enough to be in the thick-target case. To 
calculate the CR-induced 7-ray albedo of SSSBs, the Moon 
spectrum has been used as a template and scaled according to 
the size distributions and densities of bodies in different popu- 
lations, such as the MBAs, Jovian and Neptunian Trojans, and 
KBOs. The CR-induced 7-ray albedo spectrum of the Moon 
itself has been calculated by Moskalenko & Porter (2007) us- 
ing a Monte Carlo code based on the GEANT4 framework. 
Most of the Lunar 7-ray emission comes from the thin rim 
with a steep fall off with energy, essentially cutting off above 
3^ GeV. The central portion of the Moon disk has an even 
steeper spectrum with a cutoff above '^600 MeV. The model 
calculations are in excellent agreement with the EGRET ob- 
servations of the Moon (Orlando & Strong 2008; Thompson 
et al. 1997). 

For objects with sizes smaller than ~ 1 meter the CR cas- 
cade does not fully develop; this is similar to CR interactions 
with a single nucleus with the produced 7-ray flux scaling 
linearly with the total amount of material in the interaction 
column. The 7-ray spectrum in this case is harder than the 
thick-target case with a spectral slope close to that of the par- 
ent CR nuclei. The 7-ray albedo of larger bodies has a soft 
spectrum which will also contribute to the total 7-ray flux < 1 
GeV. Therefore, the total albedo spectral shape below/above 
1 GeV will depend on the relative abundances of the debris 
material and the larger rocks. 

Production of 7-rays in pp-interactions from the decay of 
neutral pions and kaons has been discussed in many papers 
(e.g., Badhwai- et al 1977; Dermer 1986a,b; Stecker 1970; 
Stephens & Badhwar 1981, and more recently Kamae et al. 
2006 and Kelner et al. 2006). We calculate the 7-ray flux 
using the method described in Moskalenko & Strong (1998) 
which is based on the work of Dermer (1986a,b). For colli- 
sions involving nuclei, the corresponding cross section is mul- 
tiplied by a factor {A^^^ + A^^^ — 1)^ where Ai and A2 are 
the beam and target nucleus atomic numbers (Dermer 1986a; 
Orth & Buffington 1976), while the energy per nucleon is the 
kinematic variable. The accuracy of this method is enough for 
the current calculations. 

Since oxygen is the most abundant element in the interstel- 
lar rock and ice we calculate the production cross section of 
7-rays for CR protons and alphas interacting with an oxygen 
target nucleus. The cross section is converted from nucleus^^ 
to g^^ units by multiplying by where Na is 

Avogadro's number The total 7-ray emissivity is obtained by 
integrating the respective production cross sections with the 
interstellar (unmodulated) incident spectrum of CR protons 




E, (MeV) 

Fig. 1. — Intensity of the IGRB as derived from the EGRET data (Strong 
et al. 2004b). Curves are shown for the thin target case for different column 
densities (top to bottom): 0.01, 0.001, 0.0001 g cm-^. 



and alphas taken from Moskalenko & Porter (2007). 

Figure 1 shows the spectrum of the IGRB derived from the 
EGRET data (Strong et al. 2004b) together with the 7-ray 
albedo of the OC debris (thin-target case) shown for different 
total column densities (top to bottom): 0.01, 0.001, 0.0001 g 
cm^^. The latter can be compared with the total gas column 
density along the line of sight to the edge of the Galactic disk 
through the Galactic center: ^0.15 g cm^^ assuming an av- 
erage 1 H-atom cm^^ and a maximum Galactocentric radius 
~20 kpc. 

The integral 7-ray emissivity in the thin-target case can be 
easily estimated: 

S{E-y > 100 MeY) ^6 X 10-^ g-^ s-\ (!) 
> 1 GeV) w5 X 10-3 g-i (2) 

and the integral intensity simply following by multiplica- 
tion with the column density, x, e.g., I{E~f > 1 GeV) = 
xS{Ey > 1 GeV) cm^^ s~^, where x is in units g cmT^. 

Any assumption of the spatial distribution of debris would 
be highly speculative. For an order of magnitude estimate we 
take the simplest case of debris concentrated in a spherical 
shell with radius d. In this case the total mass of debris is: 

Md - OAld^xM^, (3) 

where d is the distance in AU. It is enough to have ^hOM^ 
of material spread out over 47r at 10'^ AU, which corresponds 
to 10^^ g cm^^, to obtain a 7-ray albedo flux within an or- 
der of magnitude of the IGRB. This column density corre- 
sponds to ~10~'^ of the total column density along the line 
of sight through the Galactic disk. Since the d^ dependence 
is very strong, much less material is required to produce the 
same albedo 7-ray flux at a smaller effective distance and vice 
versa. 

5. DISCUSSION 

The total mass and distribution of SSSBs smaller than 
cometary nuclei are highly uncertain, as is the spatial distribu- 
tion of debris. Although, simulations indicate that the number 
of comets and the total mass considerably increase toward the 
inner edge of the OC. It is entirely possible that ~100 of 
debris is distributed in the vast space of the OC. Since the OC 
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Table 1 . Debris detection limits in different asteroid populations 



Population 


Semimajor 
axis, AU 


Total mass, 

M® 


Debris detection 
limit, 


MBAs" 


2.1-3.3 


~6 X 10-"* 


~4 X IQ-^ 


Jovian Trojans 


5.2 


~1 X 10-4 


~(1.8-3.8) X lO-*' 


Neptunian Trojans 


30 


~1 X 10-3 


~9 X 10-5 


KBOs'' 


30-50 


~0.1 


~1.3 X 10-2 


oe 


>100 


~100 


~50 



^Integral over the ecliptic longitude assuming average d ■- 
''Using d = 40 AU. 

■^Using X = 10-" g cm-2, d = 10^ AU. 



: 2 AU. 



has spherical and disk components, it is possible that the OC 
7-ray albedo is brighter around the ecliptic. 

The minimal detectable mass in debris for other systems 
can also be estimated (Table 1). For example, assuming 
that KBOs are distributed uniformly within ±10° around the 
ecliptic, the fraction of the total solid angle subtended by 
the Kuiper Belt is Cl/4n = 0.173. The detectable mass 
of KBO debris can then be calculated using eq. (3) Md = 
{n /4Tr)0. 47 cPxMs) ~ 1.3 x IQ-^Me, with x = 10"^ g 
cm"2, and d ^ 40 AU. 

Jovian and Neptunian Trojans would appear as point-like 
sources. The integral 7-ray flux from a point-Uke mass M4 of 
debris at a distance d AU is: 



F{E^ > 1 GeV) 



MdS{E^ > 1 GeV) 



0.01- 



d^ 



■ cm 



where we used the emissivity S{E^ > 



(4) 

1 GeV) correspond- 



ing to the interstellar CR spectrum calculated earlier Since 1 
GeV photons are produced by CR protons of > 10 GeV, which 



are only slightly modulated in the heliosphere, this is a good 
approximation. 

The one-year LAT integral point source sensitivity^ is 
Fiy,{E^ > 1 GeV) ~ 1 X 10-9 cm-^ s'^ This gives a 
detectability condition for one year of exposure: 



Md > 1 X 10-^Med'. 



(5) 



The proposed method is sensitive enough to detect a de- 
bris mass as small as ~0.1% - 50% of the total mass of the 
system. It is most sensitive for the MBA population due to 
its proximity, but can also provide meaningful limits on more 
distant populations, such as Jovian and Neptunian Trojans, 
KBOs, and OC. A clear signature of the described process is 
the power-law spectrum of the 7-ray albedo (>1 GeV) which 
has an index similar to that of the incident CR spectrum. The 
larger photon statistics of the projected 5-year LAT mission 
will probe even smaller debris masses for all these systems. 

The sensitivity and resolution of the LAT will allow it to 
resolve many more individual sources, such as AGN, not re- 
solved by EGRET and that contribute to current estimates 
of the IGRB. Other components of the remaining IGRB will 
therefore become more important. Understanding of the in- 
strumental backgrounds within the LAT will aUow discrimi- 
nation of the IGRB at 10% of the current level (Atwood et 
al. 2009). This will allow a meaningful estimate of the fore- 
ground from the OC and other asteroid populations to be ob- 
tained. On the other hand, correct determination of such a 
foreground will allow for more accurate determination of the 
truly extragalactic component of the diffuse emission. 

I. V. M. and T. A. R acknowledge support from a NASA 
Astronomy and Physics Research and Analysis Program 
(APRA) grant. T. A. P. acknowledges support from the US 
DOE. 
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